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used,  the  probability  of  resolving  targets  does  not  rise  above  0.9  until  the  targets  are  separated  in  range 
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RANGE  RESOLUTION  OF  TARGETS  USING  AUTOMATIC  DETECTORS 


INTRODUCTION 

There  are  two  adverse  effects  of  closely  separated  targets  on  automatic  detectors. 
First,  a target  can  lie  within  the  reference  cells  of  the  other  target,  causing  the  estimated 
threshold  to  oe  too  large  and  resulting  in  the  other  target  not  being  detected  (target  sup- 
pression). Second,  if  the  targets  are  close,  the  detections  from  both  targets  are  merged 
and  a single  target  is  reported  (targets  are  not  resolved).  The  problem  of  target  suppres- 
sion has  received  some  attention,  and  the  papers  of  Urkowitz  and  Perry  |1],  Trunk  et  al. 
[2],  and  Rickard  and  Dillard  [3]  propose  solutions  by  removing  large  returns  from  the 
reference  cells.  The  problem  of  range  resolution  of  targets  has  received  little  attention 
and  is  the  subject  of  this  report. 

The  probability  of  resolving  targets  is  not  only  a function  of  target  separation  but 
also  a function  of  the  signal  strengths.  However,  it  has  usually  been  assumed  that  if  tar- 
gets are  large  enough  to  be  detected,  they  can  be  resolved  if  they  are  separated  by  a 
pulsewidth  or,  equivalently,  lie  within  diherent  range  cells  (4|.  In  practice  this  resolution 
criterion  is  rarely  achieved.  The  basic  problem  is  that  the  shape  of  the  returned  pulse 
from  a target  is  rarely  known  in  a realistic  radar  environment.  There  are  several  reasons 
for  this.  First,  the  targets  are  not  point  scatterers  but  are  complex  scatterers.  Second, 
even  from  a point  scatterer  the  shape  of  the  returned  pulse  is  seldom  known  because  of 
uncertainties  in  the  transmitted  pulse  and  receiver  bandwidth  in  addition  to  spreading  of 
the  pulse  due  to  reflections  from  nearby  objects  such  as  buildings,  towers,  and  masts. 

The  purpose  of  this  report  is  to  investigate  and  recommend  the  type  of  video  and 
thresholding  that  should  be  used  in  automatic  detection  systems  so  that  closely  spaced 
targets  can  be  resolved.  Since  there  is  uncertainty  about  the  shape  of  the  returned  pulse, 
the  optimal  resolution  system  will  not  be  considered.  Instead  an  algonUim  that  is 
presently  being  used  in  several  systems  will  be  used  in  this  study  to  merge  multiple 
detections  from  a single  target  and  to  resolve  multiple  detections  from  different  targets. 
The  algorithm  for  a 3D  radar  is  discussed  in  the  next  section. 


RESOLUTION  ALGORITHM 

The  antenna  beam  positions  for  a typical  3D  surveillance  radar  are  shown  in  Fig.  1. 
Consequently,  depending  on  the  target  size  and  location,  the  target  can  be  detected  in 
several  azimuth-elevation  beam  positions.  Furthermore  the  target  can  be  detected  in 
adjacent  range  cells.  Thus  it  is  not  unusual  for  a large-cross-section  target  to  be  detected 
in  as  many  as  10  range-azimuth-elevation  cells.  Consequently  an  algorithm  (decision  rule) 
is  required  to  merge  these  detections  (threshold  crossings)  into  a single  centroided  detec- 
tion. 
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Fig.  1 — Antenna  beam  poaitions  for  a typical  3D 
radar.  The  circlea  represent  3-dB  contours  of  the 
antenna  beam  at  different  scanning  positions. 


The  alfforithm  selected  for  this  study  decides  whether  a new  detection  is  adjacent  to 
any  of  the  previously  determined  set  of  adjacent  detections.  If  the  new  detection  is  adja- 
cent to  any  detection  in  the  set  of  adjacent  detections,  it  is  added  to  the  set.  If  it  is  not 
adjacent  to  any  detection  in  the  set,  the  next  detection  is  considered.  Two  detections 
are  adjacent  if  two  of  their  three  parameters  (range,  azimuth,  and  elevation)  are  the  same 
and  the  other  parameter  differs  by  the  resolution  element;  range  resolution  cell  AR, 
azimuth  beamwidth  6,  or  elevation  beamwidth  7.  To  avoid  roundoff  eiron,  the  new 
detection  is  added  to  the  previous  set  of  adjacent  detections  if  any  of  three  conditions 
are  satisfied  for  any  feth  detection  in  the  set  of  K detections.  The  first  condition  is 

W?  -«*|<  1.2AR, 

\A 

lE-EfcKO.lT, 

the  second  condition  is 

\A  -A;^|<1.20, 

|E-£fc|<0.l7, 

and  the  third  condition  is 

\R-R^\<0.1AR, 

\A  -A^\<0.10, 

!£-£*!<  1.27. 
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Fig.  2 — Example  of  merging  of  adjacent 
detectiona.  The  groupings  that  are 
merged  are  indicated  by  heavy  connect- 
ing linei. 
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where  R.  A,  and  E and  R),,  Af,  and  are  the  range,  azimuth,  and  elevation  of  the 
detection  being  tested  and  the  kth  detection  in  the  set  of  adjacent  detections  respectively. 
Therefore  the  13  detections  in  Fig.  2 are  merged  into  three  centroided  detections.  The 
detections  in  the  largest  grouping  of  adjacent  detections  are 

(Aii42£i).  {R2AiE2),  («3i4,£,).  The  second 

largest  grouping  is  {R2A2E2),  (£3-4, £3),  (£3/12^3).  and  (R3A3E2).  The  final  grouping 
contains  the  single  detection  (£3i43£j). 

The  estimates  of  range  £,  azimuth  A,  and  elevation  £ for  the  K adjacent  detections 
are  ^ 
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where  Sf,  is  the  signal  power  associated  with  the  Aith  adjacent  detection. 
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Fig.  3 — Variations  of  the  basic  detector.  T is  the  threshold,  (I  either  is  the  mean  of  all  the 
reference  cells  or  is  the  smaller  (the  minimum)  of  the  means  of  the  reference  cells  at  only  the 
right  and  of  the  reference  cells  at  only  the  left,  F is  a fixed  number,  and  a is  the  variance. 


Many  types  of  detection  schemes  are  possible.  The  collection  of  detectors  considered 
is  discussed  in  the  next  section. 


DETECTION  CLASS 

The  basic  detector  uses  the  reference  cells  surrounding  the  test  cell  to  generate  a 
threshold  to  which  the  test  cell  is  compared.  Variations  of  this  basic  detector  are  shown 
in  Fig.  3.  The  set  of  detectors  uses  either  linear  or  log  video,  either  all  the  reference  cells 
or  the  half  with  the  minimum  mean  value,  and  either  a one-parameter  (mean)  threshold  or 
a two-parameter  (mean  and  variance)  threshold. 

Generally  M pulses  are  integrated  in  each  range  cell  and  beam  position,  yielding  the 
integrated  signals 

M 

' (1) 

1-1 

where  Xfj  is  the  returned  envelope-detected  signal  from  the  I'th  pulse  in  the  yth  range  cell. 
For  3D  radars,  M usually  equals  one  or  two.  Then  the  mean  value  (depending  on  whether 
the  number  of  reference  cells  Nf^  on  one  side  of  the  test  cell  are  used  or  the  total  number 
of  reference  cells  2Nd  are  used)  is  given  by 
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— I , 
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The  corresponding  mean  squares  are 


1 " 

^ “ 2Nif  ^ 
n jm  1 


z;  ■ E 
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K ,.1 


and  the  standard  deviation  is 


where  Zj  and  Zj  use  the  same  reference  cells.  The  two-parameter  threshold  is  [51 

Tj-Zj*FaOj.  (5 

regardless  of  whether  linear  or  log  video  is  being  used.  Tlie  parameter  Fg  is  used  to  set 
the  false  alarm  rate.  (The  two-parameter  threshold  for  log  video  has  dubious  meaning, 
since  the  threshold  can  be  dominated  by  the  shape  of  the  density  function  near  zero.) 
The  one-parameter  threshold  is 


for  linear  video  and 


for  log  video. 


Tj-Fg*Zj 
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Finally,  detections  are  declared  by  comparing  Zj  to  Tj,  with  a target  being  declared 
if  Zj  > Tj.  Since  it  is  difficult  to  analytically  evaluate  the  resolution  properties  of  the 
various  detectors,  simulation  techniques  were  used. 


SIMULATION 

Before  considering  the  simulation  of  the  detector,  let  us  consider  the  target  geometry. 
Two  targets  at  the  same  azimuth  and  elevation  are  separated  in  range  by  a distance  AD, 
where  1.5AR  < AD  < 3.5 AR,  with  AR  being  the  width  of  the  range-resolution  cell.  The 
signal-to-noise  ratio  for  each  target  is  20  dB  per  pulse  when  the  target  is  in  the  center  of 
a range-azimuth-elevation  cell.  Finally,  in  one  example,  a third  target  is  placed  in  the 
reference  cells,  seven  range  cells  away  from  the  first  target.  A detailed  description  of  the 
simulation  of  the  detector  subroutine  is  given  in  Appendix  A. 

In  this  simulation  the  target  was  uniformly  distributed  in  range  and  azimuth  in  the 
center  range-azimuth  cell.  The  range  cell  dimension  was  the  3-d B pulse  width,  and  adjacent 
azimuth  beam  positions  crossed  at  the  1.6-dB  point.  The  targets  were  placed  in  the  center 
of  the  elevation  beam.  Since  the  signal  energy  is  down  36  dB  in  adjacent  elevation  beam 
positions  crossing  at  the  3-dB  point,  only  one  elevation  beam  was  simulated.  Finally, 
only  one  pulse  was  used  (no  integration),  and  ten  reference  cells  were  used  on  either  side 
of  the  test  cell. 

The  detector  subroutine  declares  target  detections  by  generating  pulse-to-pulse  video 
returns,  integrating  the  returned  signal,  and  comparing  it  to  an  adaptive  threshold  which 
is  generated  from  the  surrounding  reference  cells.  In  addition  to  the  range-azimuth  cell 
in  which  the  target  is  present,  signal  energy  is  placed  in  adjacent  range  cells  (only  one  on 
either  side)  according  to  a [(sin  x}fx]^  pulse  shape  and  in  adjacent  azimuth  beams  (again 
only  one  on  either  side)  according  to  a [(sin  x)/x]^  antenna  pattern.  Thus,  in  each 
elevation  beam  position  the  target  can  be  detected  in  nine  range-azimuth  cells  (three  range 
cells  times  three  azimuth  cells). 


MONTE  CARLO  RESULTS 

The  first  case  simulated  consisted  of  the  two  targets  separated  by  1.5,  2.0,  2.5,  3.0, 
and  3.5  range  cells.  For  each  separation  ten  repetitions  were  used.  The  results  using  all 
the  reference  cells  or  using  the  half  of  the  reference  cell  with  the  smaller  mean  value  are 
given  in  Table  1.  When  the  target  separation  is  1.5AR,  almost  always  only  one  target  is 
declared  regardless  of  the  video,  reference  cells,  or  thresholding  employed.  However,  when 
all  the  reference  cells  are  used  and  the  target  separation  is  larger,  the  one-parameter 
thresholds  correctly  identify  the  presence  of  two  targets,  but  the  two-parameter  threshold 
generally  declares  no  targets.  The  targets  are  suppressed  when  two-parameter  thresholding 
is  used  because  a target  in  the  reference  cells  has  a great  effect  on  the  estimation  of  the 
noise  variance  but  has  little  effect  on  the  estimation  of  the  mean.  When  the  reference 
cells  with  the  smaller  mean  are  used  and  the  target  separation  is  2.5AR  or  greater,  all 
detectors  correctly  identify  the  presence  of  two  targets. 

The  second  case  simulated  consisted  of  two  targets  separated  by  1.5,  2.0,  2.5,  3.0, 
and  3.5  range  cells  and  a third  target  7.0  range  cells  from  the  first  target.  Table  2 shows 
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Table  1 — Number  of  Targets  Detected  in  the  Case  of  Two  Targets  Separated 
by  1.5,  2.0,  2.5,  3.0,  or  3.5  Range-Cell  Dimensions  AA 


Target 

Separa- 

tion 

No.  of  Times  a Given  No.  of  Targets  were  Detected 

Out  of  Ten  Opportunities 

Linear  Detector, 

T “ p + Fo 

Linear  Detector, 

Log  Detector, 

T-F+p  I 

Zero 

Targets 

One 

Target 

Two 

Targtte 

_aa 

Zero 

Targets 

On"  1 ’'^wo 
Target  jTaiiali 

Zero 

Targets 

One 

Target 

Two 

Three 

Targets 

Threshold  Based  on  All  Reference  Cells 

1.5AA 

2.0AA 

2.5AA 

3.0AA 

3.5AA 

mtiii 

immr 

utriiii 

imiii 

mriiii 

II 

1 

II 

1 

MlHt 

im 

1 

1 

IW 

IWilll 

wtwt 

wtwt 

wt 

1 

Wt 

Wtlll 

wtwt 

wtwt 

1 

1 

Threshold  Based  on  Minimum  of  Reference  Cells 

1.5AA 

2.0AA 

2.5AA 

3.0AA 

3.5AA 

iHtim 

iWt 

wt 

WrJHT 

iHtiW 



WtiW 

wt 

wt 

WtWf 

wtwt 

wtwt 



wtwt 

wt 



wt 

WtWf 

wfwr 

wtwt 

that  the  log  detector  using  the  reference  cells  with  the  smaller  mean  value  has  the  best 
performance  and  thus  should  be  used. 

To  determine  more  accurately  the  resolution  capability  of  the  log  detector  using  the 
minimum  of  the  reference-cell  means,  the  distance  between  the  two  targets  was  increased 
from  l.OAR  to  3.5AA  in  steps  of  O.IAA,  and  for  each  separation  40  repetitions  were 
run.  The  results  of  this  simulation  are  given  in  Fig.  4.  The  probability  of  correctly 
resolving  the  two  targets  appears  fairly  linear  from  1.7AA  to  2.6AA.  Search  techniques 
corroborated  that  the  function 


P(N  = 2IAD)  = 0, 


AD«;i.7AA, 


AD  - 1.7AA 
0.9AA  ’ 


1.7 AA  <AD<  2.6AA, 


= 1,  2.6AA<AD, 

where  P{N  = 2|AD)  is  the  probability  that  two  targets  are  detected  when  the  two  targets 
are  separated  by  distance  AD,  is  the  best  fit  whether  the  mean-square  error,  maximum 
error,  or  absolute  error  is  used  as  the  goodness-of-fit  criterion.  The  goodness-of-fit  error 
can  be  reduced  by  a factor  of  2 if  the  following  (more  complicated)  expression  is  used; 
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Table  2 — Number  of  Targets  Detected  in  the  Case  of  Two  Targets  Separated  by  1.5,  2.0, 
2.5,  3.0,  and  3.5  Range-Cell  Dimensions  A/i  and  a Third  Target  7.0  Range  Cells  From  the 
First. 


No.  of  Times  a Given  No.  of  Targets  Were  Detected 

Separa- 

Out  of  Ten  Opportunities 

tion  o^ 

Linear  Detector, 

Log  Detector, 

First  Two 

T~Fn 

T-F  + p 

Targets 

Zero 

One 

Two 

Zero 

One 

Two 

Targets 

Target 

Targets 

Targets 

Target 

TargeU 

Threshold  Based  on  All  Reference  Cells 

1.5AR 

imilii 

1 

MM 

2.0AR 

Ml 

nil 

nil 

Ml 

2.5AR 

1 

Mil 

11 

MM 

3.0AR 

ill 

III 

nil 

n 

Mill 

3.5AR 

1 

Mil 

II 

MM 

Threshold  Based  on  Minimum  of  Reference  Cells 

1.5AR 

MM 

MM 

2.0AR 

nil 

Ml 

nil 

Ml 

2.5AR 

1 

Mini 

MM 

3.0AR 

nil 

Ml 

MM 

3,5AR 

1 

III 

Mil 

MM 

P{N=2\M))  = 0,  AD<l.66AR, 

= 0.88(AD  - 1.66R)/AR,  1.66AR  < AD  < 2.16AR, 

= 1 - 0.56  exp(-  5.87(  AD  - 2.16)/AR  ] , 2.16A/?  < AD. 


CONCLUSION 

For  resolving  closely  spaced  targets,  log  video  should  be  used  and  the  threshold 
should  be  of  the  form  T = p + F where  the  mean  p is  the  smaller  of  the  two  means 
calculated  from  the  reference  cells  on  either  side  of  the  test  cell  and  F is  a fixed  number. 
This  procedure  does  not  allow  a few  large  returns  in  the  reference  cells  to  raise  the 
detection  threshold  appreciably.  Thus  this  procedure  will  not  suppress  spiky  clutter 
and  should  not  be  used  without  an  MTI  in  the  clutter  regions. 

When  the  adjacent-detection  merging  algorithm  is  used,  the  probability  of  resolving 
targets  does  not  rise  above  0.9  until  the  targets  are  separated  by  2.5  pulse  widths.  If  this 
resolution  is  not  adequate,  it  can  be  improved  by  using  other  resolution  algorithms  which 


1.5  2.5  3.5 

TARGET  SEPARATION  IN  RANGE  CELL  DIMENSIONS  AR 


Fig.  4 — Resolution  capability  of  log  detector  which  uses  reference  cells  with  the 

minimum  mean 


depend  on  assuming  a returned  pulse  shape.  If  the  pulse-shape  assumption  is  correct,  one 
should  be  able  to  resolve  targets  to  within  (approximately)  a pulse  width. 
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Appendix  A 
SUBROUTINE  DET3D 


Subroutine  DET3D  should  be  called  for  every  target  and  elevation  mode  as  long  as 
the  difference  between  the  center  of  the  elevation  beam  of  the  radar  mode  and  the  target 
elevation  angle  is  less  than  the  3-dB  elevation  beamwidth.  The  purpose  of  the  routine  is 
to  declare  target  detections  by  generating  pulse-to-pulse  video  returns,  integrating  the 
returned  signal,  and  comparing  it  to  an  sulaptive  threshold  which  is  generated  from  the 
surrounding  reference  cells.  Detections  can  be  made  in  either  adjacent  range  cells  or 
adjacent  azimuth  beam  position  in  addition  to  the  range-azimuth  cell  in  which  the  target 
is  present.  Thus,  in  theory  a radar  mode  can  report  nine  detections  of  a single  target. 

The  subroutine  initially  tests  whether  appropriate  input  parameters  are  less  than  pre- 
assigned values:  the  number  of  pulses  integrated  M is  less  than  32,  number  of  reference 
cells  Nfi  on  each  side  of  the  test  cell  is  less  than  10,  and  the  absolute  value  of  the  clutter 
correlation  coefficient  p is  less  than  1.0.  If  any  parameter  exceeds  its  limit,  an  error 
message  is  printed,  and  one  exits  from  the  subroutine. 

Next  M times  the  signal-to-clutter-plus-noise  ratio  is  compared  to  2.  If  the  value  is 
less  than  2,  one  declares  that  the  target  is  not  detected,  the  number  of  detections  of  the 
/th  target  by  the  </th  mode  [NDET3(I,  J))  is  set  to  zero,  and  one  exists  from  the  sub- 
routine. 

The  detail  simulation  is  begun  by  finding  all  the  targets  which  lie  within  the  reference 
cells  of  the  NTARth  target.  The  list  of  interfering  targets  (INF)  is  initialized  by  setting 
INF(l)  = NTAR,  and  the  corresponding  signal  power  SNREF(l)  is  set  equal  to  the  target 
power  S.  The  remaining  interfering  targets  are  found  by  calling  subroutine  RESOL  which 
calculates  the  number  of  interfering  targets  (NI),  lists  the  index  of  the  interfering  targets 
in  the  array  INF,  and  lists  the  corresponding  signal  powers  in  P. 

The  main  azimuth  do  loop,  which  generates  the  azimuth  beam  positions  on  either 
side  of  the  azimuth  beam  position  with  maximum  gain,  first  generates  the  azimuth  beam 
position  Q proceeding  the  maximum-gain  beam  position.  This  beam  position  is  given  by 

+IAZ-2),  (Al) 

where  lAZ  = 1,  6g  is  the  azimuth  angle  between  complete  elevation  scans,  and  Kg  is  the 
integer  defined  by 

Kg  = integer  [(A  - Oydg  + 0.5] , (A2) 

in  which  A is  the  target  azimuth  and  6 is  the  azimuth  angle  of  the  first  elevation  scan. 
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In  the  simulation  the  video  return  is  generated  in  only  the  reference  cells  surrounding 
the  target.  Thus,  to  save  computer  storage,  only  25  range  cells  are  saved,  and  the  target 
is  always  placed  in  the  13th  range.  Consequently  the  range  to  the  start  of  the  first  range 
cell  is 


Rg  = AR{K/fg  - 13), 


(A3) 


where  A/2  is  the  range-cell  dimension  and  Kf^g  is  the  integer  defined  by 

Kftg  = integer  [i2/A/2] , (A4) 


with  R being  the  target  range. 

The  first  step  in  generating  the  radar  video  return  is  to  generate  the  signal  (target) 
return  in  the  appropriate  range  cells.  Specifically,  the  index  of  the  flrst  range  cell  is 

(A5) 

and  the  index  of  the  last  range  cell  is 

^t  = 15  + ^R-  (A6) 

Thus,  if  Njf  = 10,  the  signal  must  be  placed  in  all  25  range  cells.  The  signal  return  for  the 
ith  pulse  and  the  yth  range  cell  (S,y)  and  an  indicator  of  signal  in  the  jth  range  cell  /(/)  are 
initially  set  to  zero  for  all  i and  ;.  Then,  the  signal  return  from  each  of  the  NI  interfering 
targets  lying  within  the  reference  cells  is  generated.  For  the  feth  target  this  is  accomplished 
by  first  generating  the  appropriate  fluctuating  amplitudes  Pik  for  the  ith  pulse,  if  the  valid 
fluctuating  amplitudes  have  not  been  calculated  previously:  IKEY(NTAR)  = 1 indicates 
that  fik  has  been  calculated  previously  for  Swerling  cases  1 and  3 (scan-to-scan  fluctua- 
tions) for  either  a previous  azimuth  or  elevation  beam  position.  If  Fik  needs  to  be 
calculated  (either  the  first  time  for  Swerling  cases  1 and  3 or  every  time  for  Swerling 
cases  2 and  4),  is  given  for  the  appropriate  Swerling  case  (NSW)  by  the  following,  in 
which  the  Uj  are  independent  random  variables  uniformly  distributed  between  0 and  1: 

• For  NSW  = 0 (nonfluctuating  target) 

= (=1 M;  (A7) 

• For  NSW  = 1 (scan-to-scan  fluctuations,  Rayleigh  density) 

= -log(f/j),  i=*l M;  (AS) 

• For  NSW  = 2 (pulse-to-pulse  fluctuations,  Rayleigh  density) 

Fik  - - log  (C/,),  i - 1, ...,  Af;  (A9) 

• For  NSW  = 3 (scan-to-scan  fluctuations,  chi-square  density) 

Fi*  =-0.5[log((/i)  + log((/2)l,  i*l M; 


(AlO) 
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• For  NSW  > 4 (pulse-to-pulse  fluctuations,  chi-square  density) 

F,*  = - O.SIlog  (I/,)  + log  , » - 1 Af.  (All) 

The  signal  is  next  placed  in  the  appropriate  range  cell  and  the  adjacent  range  cells 
by  reducing  the  returned  signal  by  a [Uin  x)lx\^  pulse  shape  and  a [(sin  x)lx\* 
antenna  pattern.  The  return  signal  from  the  kth  target  is  centered  in  range  cell  K/f, 

Kji  = integer  [(F*  -RgUAR],  (A12) 

where  Rh  is  the  range  of  the  feth  target.  The  signal-return  reduction  (F)  in  the  adjacent 
range  cell  K-p, 


Kp=Kn*I,  /--I,  Q,  1 (A13) 

(because  of  the  (sin  x)lx  pulse  shape),  is  given  by 

F=  [(sinFd)/Fdl2,  (A14) 

where 

Fa  = 2.7832(fife  - Rt)I^-,  (A15) 

and 

Rp  = (Kp  + 0.5)/!iR  + Rg.  (A16) 

At  this  time  the  indicator  HKp)  is  set  equal  to  1.  Similarly  the  signal  reduction  (G)  due 
to  the  (sin  x)/x  antenna  is  given  by 

G = [(sinGd)/Grfl4,  (A17) 

where 

Grf  = 2.7832(Afe -0)/03dB.  (A18) 

in  which  03()b  is  the  antenna  3-dB  azimuth  beamwidth  and  is  the  azimuth  of  the  feth 
target.  Finely  the  signal  (normalized  by  the  clutter  power  C and  noise  power  N)  due  to 
the  feth  target  in  the  Kp  range  cell  (for  the  azimuth  ^am  position  specified  by  (Al))  is 

S.  ^new)  = S.  ^old)  + (G){F)(F,.fe  )Pfe  /(C  + N).  ( A19) 

The  calculation  indicated  by  (A13)  through  (A19)  is  first  repeated  for  the  adjacent 
range  cells  indicated  in  (A13).  Then  the  calculation  indicated  by  (A7)  throu^  (A19)  is 
repeated  for  all  NI  targets  in  the  reference  cells.  Thus,  at  the  end  of  all  the  above  repeti- 
tions, S^  j is  the  signal  power  in  the  ith  pulse  and  ;th  range  cell  due  to  all  the  targets  in 
the  reference  cells. 
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Next  Rayleigh  noise  (and  possibly  correlated  clutter)  is  added  to  the  signal  to  pro- 
duce the  total  video  return  X^.  The  video  return  (because  of  computer  speed  considera- 
tions) is  generated  for  three  dutinct  cases:  clutter  insignificant  and  no  signal  present  in 
the  range  cell,  clutter  insignificant  and  signal  present  in  the  range  cell,  and  clutter  signifi- 
cant. The  significance  of  clutter  is  indicated  by  (C)(Af)  being  greater  than  N,  and  signal 
present  in  the  7th  cell  is  indicated  by  /O)  “ 1.  Thus  the  ith  return  in  the  ;th  ceil  (Xy)  is 
given  by  the  following: 

• When  clutter  is  insignificant  ((C)(Af)  < N)  and  no  signal  is  present  in  the  yth  cell 

m = 0), 

Xij  = (-  2 log  C/,)1'2,  , = 1, Af,  (A20) 

where  the  are  independent  uniform  random  numbers  (0.1),  different  for  each  j; 

• When  clutter  is  insignificant  ((C)(Af)  < N)  and  signal  is  present  in  the  yth  cell 

m = 1), 

Xij  = {[a,  cos  di  + y/2  (Si,)i/2]2  + (a,-  sin  0,)2}i'*,  (A21) 

where 

a,  = [-21ogt;;i^/^  (A22) 

and 

di  = 2nUi^M;  (A23) 

• When  clutter  is  significant  ((C)(Af)  > N), 

Xij  = {[CX,  + a,  cos  0,  + (2S^)1/2]2  + ^ (A24) 

where  the  Rayleigh  noise  (a,-  cos  0,-  and  a,-  sin  0,)  is  given  by 

a,.  = N'(-21ogU,>M)^'^  (A25) 

<t>i=2xUi^3M>  (A26) 

and 

S’  = [NI{C  + N)]^I^,  (A27) 

the  clutter  is  given  by 

CX,  = pCX,_i  + (1  - p2)i/2  a'i  sin  6-  (A28) 


and 
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in  which 

a;-C'(- 2 log  I/,)  *'2  (A30) 

and 

0;-27rt/,*2M.  (A31) 

with 

C'  = (C/(C  + N)1‘'2,  (A32) 

and  the  initial  clutter  values  are 

CXq  = C'ttQ  cos  00  (A33) 

and 

CYq  = C'uq  sin  Oq.  (A34) 

Next  a decision  is  made  whether  to  use  linear  or  log  video.  Then  the  M pulses  are 
integrated  in  each  range  cell,  yielding  the  values 


M 

Z;  ■ E • 
1-1 


(A35) 


which  equation  was  given  in  the  main  text  as  (1).  The  detection  threshold  Tj  for  the  y'th 
range  cell  0 = 12,  13,  and  14)  uses  either  all  the  reference  cells,  the  minimum  of  the 
reference  cells,  or  the  maximum  of  the  reference  cells.  Furthermore,  the  threshold  may 
be  either  a one-parameter  (mean)  threshold  or  a two-parameter  (mean  and  variance) 
threshold.  Thus  the  mean  is  given  by  (2),  (3),  or  (4).  The  corresponding  mean  squares, 
the  standard  deviation,  the  two-parameter  threshold,  and  the  one-parameter  thresholds  are 
given  by  (5)  through  (11). 

Finally,  detections  are  declared  by  comparing  Z12,  ^13.  and  to  T12,  Tj3,  and 
Ti4  respectively.  If  Zi  is  greater  than  Tj,  a detection  is  declared  in  the  ;'th  range  cell,  the 
counter  for  number  of  detections  (parameter  labeled  II)  is  increased  by  one,  and  the 
following  detection  parameters  are  saved: 


II 


Rs  * U + 0.5)AR 
9 

Zj(C  + N) 
Time 


the  number  of  detections, 

the  range  of  the  detection, 

the  azimuth  of  the  detection, 

the  signal  amplitude  of  the  detection, 

the  time  of  the  detection. 


Also,  if  any  target  lies  within  three  range  cells  and  2.4  azimuth  beamwidths  of  the 
NTARth  target,  MER3(NTAR)  is  set  to  -1,  noting  this  target  interfering  condition. 
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After  the  detection  tests  have  been  performed  for  the  initial  azimuth  beam  position 
(lAZ  = 1),  all  calculations  are  repeated  for  the  other  two  beam  positions.  If  there  are 
target  interfering  conditions  and  if  the  target  has  not  been  detected  on  this  mode  or  pre- 
vious modes  of  this  radar  scan,  MER3(NTAR)  is  set  to  -1.  Then  one  returns  to  the  main 
program. 
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